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Abstract 

In this paper, we will develop an iterative procedure to determine the detailed asymptotic 
behaviour of solutions of a certain class of nonlinear vector differential equations which approach 
a nonlinear sink as time tends to infinity. This procedure is indifferent to resonance in the 
eigenvalues. Moreover, we will address the writing of one component of a solution in terms of 
the other in the case of a planar system. Examples will be given, notably the Michaelis-Menten 
mechanism of enzyme kinetics. 
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1 Introduction 
1.1 The Problem 

Suppose we are presented with a vector ordinary differential equation with time as the independent 
variable. Suppose also that the origin is an asymptotically stable equilibrium point, that is, all 
solutions starting sufficiently close to the origin will approach the origin as time tends to infinity. 
It is very desirable to obtain the asymptotic behaviour of solutions as time tends to infinity with 
sufficient detail when starting close to the origin. In this paper, we will develop an iterative 
procedure, which is indifferent to resonance in the eigenvalues, to construct asymptotic expansions 
of solutions with any desired accuracy. This detailed asymptotic behaviour, for example, will 
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enable us to construct phase portraits with more detail than the linearization since the iterative 
procedure yields higher-order (nonlinear) approximations. Moreover, the detailed behaviour will 
have implications for concavity results, parameter estimation and reduction methods in chemical 
kinetics, and it will show the limitations of using traditional power series. 

Let n € N be the dimension and consider the constant matrix A € M"^" and the vector field 
b : M" — >■ M". In this paper, we will be considering the general initial value problem 

i = Ax + b(x) , x(0) = xo, (1) 

where t is time and ' = d/dt, when A and b satisfy certain properties. Note that ([1]) is equivalent 
to the integral equation 

x(t) = 6*^X0 + / e(*-^)^b(x(s)) ds. (2) 
Jo 

The matrix A is assumed to be Hurwitz with eigenvalues {Aj}"^j^ having respective real parts 
{/^illLi satisfying the ordering 

/in < fJ-n-i < ■ ■ ■ < fJ-i < 0. (3) 
For the vector field b, we assume there are constants 5,ki,k2 > such that b G C^{Bs, K") with 

||b(x)|| < A;i |lx||" and ||Db(x) || < fca ||xf for all x G (4) 

for some a > 1 and /3 > 0, where 

'56i(x) " 



Y^x] for X G M", 05 := {x G : ||x|| < (^} , and Db(x) 

\ i=i 



/i,i=i 



are, respectively, the Euclidean norm (chosen arbitrarily), the open ball of radius 6 centred at the 
origin, and the Jacobian matrix. 

Remark 1. The matrix A, the vector field b, the constants and {/ii}"=i, and the constants 

a and (3 retain their meaning throughout the paper. However, the constants 6, ki, and ^2 do not. 
We will be dealing with many combinations of estimates which are introduced by saying something 
like "there are some constants 5,k > such that..." and, typically, we take d > small enough so 
that multiple estimates apply. 

We will denote by (^j(xo) the flow of ([1]). Where necessary, we will write <^(t,xo) with compo- 
nents {(/)i(t, xq)}"^^. Furthermore, define the important ratio 

K — — 

Observe that, by virtue of ([3]), k > 1. We will need to consider two distinct cases. 

Definition 2. If 1 < k < a, the eigenvalues of the matrix A are closely-spaced relative to the vector 
field b. If K > a, the eigenvalues of A are widely-spaced relative to b. 
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The goal of this paper is to find the detailed asymptotic behaviour of the flow (^^(xo) as t — )• oo. 
The techniques we will develop are indifferent to resonance in the eigenvalues of A. However, it is 
important to be clear what resonance is and when it occurs. There is resonance in the eigenvalues 
if there exist {rui}^^^ C Nq and j G {1, . . . , n} with Y17=i ^ 2 and Xj = Yl^=i "T-i-^i- The order 
of the resonance is Y27=i easy to verify that resonance can only occur if k > 2. 

In the remainder of this section, we provide some background information and establish impor- 
tant estimates. In ^j2l we focus on the case where the eigenvalues are closely-spaced relative to the 
nonlinear part. To approximate the solution 0t(xo), we will construct iterates {^^'"•'(t, Yo)}^^^- 
The difference (in norm) of the first iterate X^^^'tiYo) the flow ^t(xo) is, for any a > 0, 
0^go(^i+(T)t iIxqII") as t — 7- oo and ||xoj| 0. Moreover, each iteration increases the closeness by 



agonalizable, and the eigenvalues are widely-spaced relative to the nonlinear part. We construct 
iterates, similar to the closely-spaced case, to approximate the flow. We also develop a result on the 
expression of the second component (j)2{t,XQ) of the flow in terms of the first component 0i(i,xo) 
when there is resonance and a quadratic nonlinearity. The general n-dimensional case is presented 
as well. In ^ we apply our results to the Michaelis-Menten mechanism of an enzyme-substrate 
reaction to generalize a result in the authors' [6]. Finally, in ^ we state some open questions. 



1.2 Background 

The two main contributions of Poincare relevant to this paper are normal forms, which originated 
in his Ph.D. thesis, and a result often referred to simply as Poincare's Theorem, which he proved 
in 1879. These are contained in the first volume (of eleven) of the Oeuvres of Poincare |22j . 

Normal-form theory involves the substitution of an analytic, near-identity transformation to 
simplify the nonlinear part of an ordinary differential equation. Resonance in the eigenvalues of A 
limits how simple the nonlinear part can be made. See, for example, [H [U |2Tl I24j . 

Consider the two systems x = Ax + b(x) and y = Ay, where b is analytic in a neighbourhood 
of the origin and satisfies ||b(x)|| = 0^||x||^^ as ||x|| — t- 0. If A has non-resonant eigenvalues and 
the convex hull of the eigenvalues does not contain zero, then Poincare's Theorem says that there 
is a near- identity, invertible transformation h with both h and analytic at the origin such 
that h((/)((xo)) = e*^h(xo). That is, the flows of the two systems are analytically conjugate in a 
neighbourhood of the origin. This helps greatly in determining qualitative properties of solutions 
of more complicated systems. Furthermore, if we write h~^(y) = y + r(y), then the relationship 
0t(xo) = e*'^h(xo) + r(e*'^h(xo)) can be used to extract asymptotic expansions. 

The Hartman-Grobman Theorem is a conjugacy result similar to Poincare's Theorem. See, 
for example, [H [71 [HI [12l [131 [13 • See also [3] which addresses the issue of just how small the 
neighbourhoods of the origin must be. Conjugacy results such as the Hartman-Grobman Theorem 
(and its variations) and Poincare's Theorem are distinguished by the specifics of the system: The 





In ^ we focus on the special case where n = 2, A is di- 
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smoothness of b, the smoothness of h, whether or not the eigenvalues of A are resonant, and 
whether or not the origin is a hyperbolic fixed point. 

The standard methods have well-known limitations. Many results have to exclude resonance and 
those that do not are often limited in their practicality. Leading-order behaviour of solutions can 
be difficult to obtain and more detailed behaviour can be more difficult still. Finally, the standard 
methods often give little or no detail on the relationship between components of solutions. 

In this paper, we will frequently encounter a certain class of functions defined by improper 
integrals. Let 5 > be a constant and define the set 

fts := { (i, x) : t > 0, X G M", ||x|| <6}. 

Let f € C{fls,W^) and suppose there are k,p > such that ||f(t,x)|| < ke^f^ for every (t,x) G Qs- 
Then, g(t,x) := f(s,x) ds, defined for (t,x) G $7^, satisfies ||g(t,x)|| < k/p for all (t,x) G ^s- It 
follows from the Weierstrass M-Test that the integral converges uniformly on fls and g G C (fi^, R"). 
See, for example. Theorems 14-19 and 14-22 of [1] and §6.5 of |19j . 

1.3 Important Exponential Estimates 

Claim 3. For A := diag{Xi, . . . , Xn), = ^'^'^ and ||e^*^|| = e^^"* = e^^^^* for every t > 0. 

Proof: The proof is straight-forward and omitted. □ 

Lemma 4. Consider the matrix A. For every a > 0, there exist ki,k2 > such that 

||e*^|| < ki e(^i+'^)* and ||e-*^|| < fca e(-^'"+'^)* = yfcs e(^'^^i+'")* for all t > 0. (5) 
Moreover, if A is diagonalizable, then we may take a = in 

Proof: The proofs of the estimates when a > are standard. Assume A is diagonalizable. That 
is, there is an invertible matrix P such that A = PAP~^, where A := diag(Ai, . . . , A^). Thus, 
gtA _ Pe*-^p^i. The estimates when a = then follow from Claim [3l □ 

Theorem 5. Consider the flow (f)t{xQ) for ([TJ. For any a G (0, —pi), there exist 6,k > such that 

||0t(a^))|| < A;e(^'i+'")* lla^ll for all {t,xo)£ns. (6) 

Proof: We will assume that ||xo|| > for if this were not the case the result would be trivial. This 
proof is more or less standard but we are careful to include the initial condition in the estimate 
which is not standard. See, for example, pages 314 and 315 of [9j. In [9j, they used Gronwall's 
Inequality at a key step (they had a slight difference in assumptions) but we will instead use a 
modification. Such modifications are examples of Bihari's Inequality. See, for example, [2]. 
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It follows from the integral equation ([2]) and the estimates ([!]) and ([5]) that 

||<^t(xo)|| < A:ie('^^+-)* ||xo|| + A^s T e(^i+-)(*-«) ||<^,(xo)r ds 

Jo 



(7) 



for some 5i,ki,k2 > 0, provided ||(?!)t(xo) |1 < di for all t>0. Later, we will show this condition is 
satisfied for ||xo|| sufficiently small. Manipulate ([7]) to obtain u{t) < v{t), where 

u{t) := e-(^'i+'^)* \\M^o)\\ > and v{t) := h ||xo!| + ^2 / e("-^)(^i+'^)"u(s)" ds > 0. 

Jo 

Now, v{t) = A;2e(""^)(^'i+'^)*'u(t)" > 0. Since u{t) < v{t) and v{t) > 0, we have 



< k2 e("-i)(^i+'^)«, v{0) = ki llxoll for s > 0. 



(8) 



If we integrate ([8]) with respect to s from to t and recall that a > 1 and /ii + cr < 0, we obtain 



1 



1 



1 — a 
Rearranging 



v{t)'^-'^ llxoir"^ 



< 



ko 



{a - 1) {fii +a) 



< 



ko 



(a - 1) (/ii + cr) 



1 



> 



l^ol 



IQ-l 



_! k2 ||xo 



la-l 



/il + fJ 



(9) 



Now, we want the expression in the brackets to be strictly positive, which yields the condition 



Take 



<5 : = 



P_ 

k^ 



ko 



|xo|| < 



and k :- 



ki 



1 k2S^ 
k^ ^ fii+a 



' " = A:i (1 - /)°-^) > 0, 



where p G (0, 1) is arbitrary, and assume ||xo|| < S. Note k is independent of ||xo||. From 

1 1 



> 



|X0 



Rearranging and recalling that u{t) < v{t), we have u{t) < k ||xo||. By definition of u{t), we have 
\\M^o)\\ < A;e(^i+'^)* llxoll for all t > 0. 

All that remains is to show that the condition ||0t(xo) || < 5i for alH > is satisfied for a partic- 
ular choice of p. Using our expressions for k and 5, we see that we can choose p £ (0, 1) sufficiently 
close to zero to ensure k5 < 6i. Hence, if (t,xo) E fls then \\4>ti'^o)\\ < ke^^^^'^'^^^ ||xo|| < k5 < 5i. □ 

Remark 6. A glance at the proof of the previous theorem shows that a is only necessary in the 
basic decay rate of (pti'^o) if is necessary in the decay rate of e*-^. Indeed, if A is diagonalizable 
then we can take o" = in (161). 
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2 Closely-Spaced Eigenvalues 

Consider once again the system ([T]). In this section, we wih assume that the eigenvalues {Aj}"^^ 
of A are closely-spaced relative to the nonlinear part b(x). That is, in addition to ([3]) and @ 
holding, we are assuming 

K<a. (10) 
Remark 7. Throughout this section, if A is diagonalizable we can take u = 0. 

Remark 8. Consider the scalar case n = 1 with x = ax + b{x), where a < 0, so that A = (a) and 
K = 1 < a. Trivially, the results of this section apply and we can take £7 = for all estimates. 

2.1 The Transformation 

Claim 9. For any a > there are 5,k > such that 

||e-*^6(0t(a^))|| < fce[("-'^)'^i+'^l*||so!r for all {t,Xo)ens. (11) 

Proof: It follows from (gD, dS]), and ([6]). □ 
Recall the integral equation When ||xo|| is sufficiently small, (jlOp and (jlip allow us to write 

M^o) = e*^ xo + / e-^^b(0,(xo)) ds - / e(*-^)^b((/.,(xo)) ds. (12) 
Jo J Jt 

See, for example, the proof of Theorem 4.1 (Stable Manifold Theorem) of Chapter 13 of [S] which 

involves the trick of "flipping an integral" in an integral equation. Now, there is a 5 > (the one 

in pip will work) such that we can define the function 

/■oo 

^(xo) := XO + / e-^^b((/)3(xo)) ds, xq G Bs- (13) 
Jo 

Using ()12p . we have thus proven the following. 
Claim 10. There is a 5 > such that 

e^'-'^'^biMxo)) ds for all {t,xo)ens. (14) 

Claim 11. There are 6,k > such that -0 € C{Bs,W) and 

|j'0(a^) — a^ll < /c lla^ll" for all xq e Bs- (15) 
Moreover, xjj is a near-identity transformation. 



Proof: It follows from (|lip . (|13p . the continuity of the integrand, and the fact that a > 1. □ 

Remark 12. The near-identity transformation i/> has many interesting properties which we will 
cover in a later paper. For example, ■j/'(xo) is the optimal (as t oo) initial condition for the 
linearized system y = Ay. Moreover, i/r satisfies the conjugacy condition ■0(0j(xo)) = e*'^'j/'(xo) 
and, under modest assumptions, is as smooth as the nonlinear part b. 
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2.2 Iterates 

2.2.1 Definition 

Taking inspiration from the integral equation (jl4p . we will construct iterates {x'-™'*(t,yo)}^_^ to 
approximate the flow <^f(xo). Assuming ||yo|| is sufficiently small, take as the first iterate 

x(^)(t,yo):=e*Vo (16a) 
and define the remainder recursively via 

X^"'^'Ht,yo) ■■= e*^yo - e(*-«)^b(x(")(s,yo)) ds (m G N) . (16b) 

To connect the iterates with the flow, we take yg := 'j/'(xo), where ||xo|| is sufficiently small. With 
this choice of yg, the first iterate is the best linear approximation (as t — )• oo) of (^f(xo). 

Remark 13. The computation of the iterates {x^™^(t, yo)}^^^ does not require us to know yg, 
which can be treated as a parameter, in terms of the initial condition xq. The resulting iterates, 
in conjunction with Theorem 1151b elow . tell us the form of the asymptotic expansion for the actual 
solution 0j (xq). However, it is possible to obtain an approximation of any desired order — sometimes 
even the exact expression — for yQ in terms of xq. See ^2.3.3i 

2.2.2 Existence, Decay Rate, and Closeness to the Flow^ of the Iterates 

Proposition 14. Consider the transformation defined by (jl3p and the iterates defined by ()16p . 
Let a > be small enough so that (a + 1) o" < {k — a) fii. Then, there is a 5 > (independent of 
m) such that if \\xq\\ < 6 then x^^\'ti'4'{'^)) exists for each m G N. Moreover, there is a k > 
(independent of m) such that 

X^"'\t,ip{xo)) < k e^f^'^"^^ \\xo\\ for all {t, xq) e Cls, m e n . (17) 

Proof: Let yg := i/'(xo). We need to be careful to find 5,k > which work for each m. 

o Let 5i,ii > be such that ||b(x)|| < ii ||x||" for all x G Bsj- See Equation ([3]). 

o Let ^2 > be such that ||(^f(xo)|| < 6i for all (t,xo) G il^j- See Equation ([6]). It follows that 

if ||xo|| < 52 then yg exists, 
o Let £3 > be such that ||yo|| < ^3 ||xo|| for all xq G Bg^. See Equation ([15]). 
o Let £4 > be such that ||e*^|| < 4e(''i+'')* for all t > 0. See Equation (P. 
o Let 4 > be such that ||e-*^|| < 4 e("'"^i+'")* for all t > 0. See Equation (P. 



There is no harm in taking 5i small enough so that 
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Note that the restriction on a impHes 

/xi < //I + (T < and {k — a) — {a + 1) a > 0. 

Take 

k := l^U 



(18) 



1 



a-l 



-, -1 



> ^3^4 and (5 € ( 0, min <! (^2 



K—a) ni — {a+l) a 

The upper bound imposed on ensures k is defined and greater than ^3^4. Observe 



(K — a) /ii — (q+1) (t 



k. 



(19) 



The proof of the result will be by induction on m. Assume that {t, xq) € $^5. 

Consider first the base case m = 1. Since yg exists, so does X^^Ht^yo) = ^^'^Yo- Furthermore, 

X^^Ht,yo) < ^4 6^'^^+'^)* llyoll < 44e(^i+'^)* < fce('^i+-)* |lxo|| . 

Thus, the result is true for m = 1. 

Now, assume the result is true for m € N. That is, X^™'^(t,yo) exists for xq € Bs and satisfies 
||x*-'"'*(i, yo)|| < A;e(^i+°")* ||xo||. Since ||xo|| < 6, we see ||x (™)(i,yo)|| < <5i and thus x("^+^)(t, yo) 
exists. Using ()16bp . the induction hypothesis, and the estimates given at the beginning of the proof. 



x^"'^^Ht,yo) < ^34 6^'^^+^^* llxoll +£i4fe"e(''^i-'^)* 



J{Q!-k)Mi+{Q!+1)o-]s 



|xo| 



(K—a) /xi — (a+l) (T 
Since < e("~i)('^i+'^)* < 1 and ||xo|| <6 < 6i/k, we have 



(K—a) fii — {a+l) a 



Using (fT9]l . ||x^™'^^^(i,yo)|| < fce^^'^^'^^* ||xo||- Hence, the result is true for m + 1. By induction, 
the result is true for each m € N. □ 



Theorem 15. Consider the transformation defined by (jl3p and the iterates defined by (jl6p . 
-Lei o" > 6e small enough so that (a + 1) a < (k — a) fii. Then, there are constants 6 > and 
{km}m=i *^ ^+ such that, for all {t, xq) G ils and m G N, 



|a+{m-l)/3 



(20) 



Proof: Let yg := ■j/'(xo). We need to be careful to find 5 > which works for each m € N. 

o Let 5i,ii > be such that ||b(x)|| < £1 ||x|[" for all x G Bs^- See Equation (jl]). 
o Let 62,(2 > be such that ||Db(x)|| < £2 \\^\f for all x € Bs2- See Equation 1^. 
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o Let 63,13 > be such that ||</>t(xo)|| < ^36(^1+'^)* ||xo|| for all (t,xo) G Qss- See Equation ([6]). 
o Let 6i,U > be such that llx^'"^*^ yo)|| < 4e('^i+'^)* ||xo|| for all (t,xo) G fls^ and m G N. 
See Equation pT]) . 

o Let 4 > be such that ||e-*^|| < 4 e(-'"^i+^)* for all t > 0. See Equation 
o If ||xo|| < min {53,54}, then any point on the line segment connecting <^f(xo) and ^'■'"''(i, Yq) 
is bounded in norm by {£3 + £4) e('^i+°")* ||xo|| for all t > 0. 

The proof of the theorem will be by induction on m. Take 



5 := min 



Si 61 



£4 £3+ 1 



■,53,(54 



Hence assume (t,xo) G i^s so all necessary estimates apply. 

Consider first the base case, m = 1. Using (|14|) and (|16ap along with the estimates at the 
beginning of the proof, 



poo 

-I 









|jb(0,(xo))|| ds< 



(k— a) /ii — (a+1) fj 



Using (|18p . which is valid in this proof as well, ki > 0. Thus, the theorem is true for m = 1. 

Now, assume the theorem is true for fixed m G N. First, for any s > consider the expression 



b((/.,(xo)) - b(x(™)(s,yo))|| < i2 (^3 + 4)^e/^(^^+-)^ ||0,(xo) - x('")(s,yo) 



|xo| 



where we applied the Mean Value Theorem and estimates given at the beginning of the proof. By 
the induction hypothesis, we know 



b(0,(xo)) - b(x('")(5,yo))|| < £2 {£3 + £4fe^'^^'+'^> {A:™e["+( 



m-l)/3][AJi+o-]j 



|xo 



|a+(m-l)/3l 



ixol 



= k„,£2 {£3 + £if e("+™^)(^i+'^)^ 11x0 ir+"''' . 
It follows from the integral equation p4p for 0t(xo) and the definition ()16bp for X^^^^^^^iYo) that 



(m+l) 



(i,yo) 



/>oo 













b((/,,(xo))-bfx(-)(s,yo) 



< kra£2 (£3 + £if £5 \ / e('^'^i-'^)(*~")e(°+'"^)(^i+'")' 



|xo 



km£2 {£z + UY^^ 



(k — a) /ii — (a + 1) 0" — m/3 (;Ui + u) 



g(a+m/3)(^Ji+(T)t |j^^||a+m^ 



Using psp . fem+i > 0. Thus, the theorem is true for m + l and, by induction, for all m G N. □ 
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Remark 16. Consider the constants {km}m=i ^+ from Theorem 1 1 5 1 and its proof. It is apparent 
from their construction there is some constant r > 0, independent of m, such that fcm+i ^ (r/m) km 
for all m S N and hence k^+i < (r™/m!) ki for ah m € N. Consequently, km — > as m — )• oo. 

Corollary 17. Consider the transformation defined by ()13p and the iterates defined by (|16p . There 
exists 6 > such that lim i/'(a^)) = cf)t{xo) uniformly in Qs- 

Proof: It follows from Theorem [15] and Remark [T6j □ 

2.3 Approximations 

2.3.1 Truncating the Iterates 

Theorem 1151 guarantees that X*-'"H^>yo) correctly approximates (pti^o) up to a certain order when 
Yq = tp{xo). However, yo) may contain higher-order terms. Fortunately, as is evident from 

the proof of the theorem, we can discard the terms of ^'•'"^(i, yo) which are not guaranteed to be 
correct by the theorem. If we use this simplified iterate in the calculation of x'-'^'^^H^i Yo) '^i^l 
not sacrifice accuracy. This potentially can simplify the calculation of iterates greatly. 

2.3.2 Truncating the Nonlinear Part 

The integration in the recursive definition ()16bp of the iterates may be difficult if not impossible. 
However, it is possible to replace b(x) with a Taylor polynomial of an appropriate order and not 
sacrifice accuracy. Assume that b € C^{Bs, M"), where 5 > is small. By Taylor's Theorem, we can 
write b(x) ~ Y1'S=£ as ||x|| for some index i G {2, 3, . . .}, where each bj(x) has components 
which are homogeneous polynomials of degree i. Note that we can take a = £ and f3 = i — 1. Define 
b^''"''(x) := Yl^e^ii^) each m £ {£,£ + 1, . . .}, which is the Taylor polynomial of b(x) of order 
m. Assuming ||xo|| is sufficiently small and yg = we can define X*-™''''^^^; Yo) ™ exactly 

the same way as before except with b(x) replaced by b^["^'^^l[^~^]-'(x). The resulting closeness of 
x(™')(t,yg) to <^f(xo) is no different than that stated in Theorem [T5l 

2.3.3 Approximating the Transformation 

Consider the pivotal transformation ip, defined in l\VS\] . The iterates p6|) can be computed in terms 
of Yq without knowing ■j/'(xo) explicitly. Moreover, the iterates give the form of the asymptotic 
expansion for 0((xo) with error estimates courtesy of Theorem [T5j However, we can approximate 
if necessary with an iterative method. In a later paper, we present an alternative method 
for finding or approximating ■j/'(xo). 

We will construct a sequence of approximations {iP^"^\'Xq)^'^^-^ for ■j/'(xo). Assuming ||xo|| is 
sufficiently small, Equations (fT3]) and (fT5|) along with Theorem [15] suggest we take 

V'(i)(xo):=xo and V^"^+^^(xo) := xq + e-«^b(x(") (.s, v(™)(xo))) (m G N) . (21) 
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We can state how close tp^'^\'Xo) is to i/>(xo). Claim [THl required to prove Proposition [T9l can be 
combined with Theorem 1151 and Proposition [19] to state how close x^™"^ (t, i/>*'™^(xo)) is to (^^(xo). 

Claim 18. Consider the transformation defined by (|13|) and the iterates defined by (|16|) . For any 
sufficiently small a > 0, there exist 5,k > such that, for all xq, j/q E Bs, t >0, and m gN, 

:(")(t,^(a=o)) - X^"'Ht,yo)\\ < A;e(''i+'^)* mxo) - y^W . 
Proof: The result follows intuitively from (|16p and the proof is omitted in the interest of space. □ 



Proposition 19. Consider the transformation defined by (jl3p and the approximations defined by 
(pTj) . There are constants 6 > and {A;m}m=i ^+ swc/i that 



^{xo)-'iP^^Hxo) <k,n\\xo\r'"'"'' for all xo e Bg 



Q+(m-l)/3 



and lim i]j^"^\xq) = ijjixQ) uniformly in Bs- 



Proof: The proof is omitted in the interest of space. It involves induction, the Triangle Inequality, 
the Mean Value Theorem, Theorem 1151 a-nd Claim [181 D 

2.4 Examples 

2.4.1 Star Node with Quadratic Nonlinearity 

Take 

A := h and b(x) := f ^^"i + + ^^^^i) , 

\0 a J \b2nX1 + h2\2X\X2 + 62222:2 / 

where a < and {bijk\\ji.^Y ^ ^ constants. Assume that at least one is non-zero and 
consider the system x = Ax + b(x). The origin, as we see, is a star node. Now, e*-^ = e°*I, 
/ii = a = (J.2, and k = 1. Also, we can take a = 2 and /3 = 1. Thus, k < a and the eigenvalues are 
closely-spaced relative to the nonlinear part. 

Define yg := ■'/'(xq) and x(t) := 0t(xo), where ||xo|| is sufficiently small. The first iterate, 
which gives us the linearization, is X^^\tiYo) = 6"*yo- Theorem [T5] tells us xi{t) ~ 6*^*7/01 and 
X2{t) ~ e'^*yo2 as t ^ 00. Provided yoi 7^ 0, X2{t) ^ (2/02/2/01) 2:1 (t) as t ^ 00. The Hartman- 
Grobman Theorem yields basically the same conclusion. We will use the second iterate to obtain a 
more detailed relationship between xi{t) and X2{t). Using (]16bp and the fact b(e°*yQ) = e^°*b(yQ), 

x(2)(t,y^) = e'^Vo + e2'^*e where ^= b(yo) . (22) 
o Case 1: yoi = 2/02 = 0. If yg = 0, X^^^^^yo) = for each m € N. By Corollary [JTl x(t) = 0. 
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o Case 2: yoi 7^ 0. It follows from (j22p and Theorem [T5l that, as t — > oo, 

xi{t) = e'^'yoi + e^'^^Ci + O {e^"') and X2 (t ) = e°*yo2 + e^^^^^ + O (e^''*) . (23) 
Using ([23]), xi = 0(e"*) as t oo and e"* = 0{xi) as xi ^ 0. Now, the first equation of ([23]) 



can be written, as t — )• oo 
at I ^1 



yoi 



1 



yoi. 



Vol) WmJ 



where we used the fact [l + M + 0(ti2)] ^ = l-u + 0{ v?^ as ti — 7> 0. Back-substituting e°*, 

^" Ki^) - (I) ^ ''(^'"*' ^ (;;^) - (I) ^' ^ 



.^01/ vyoiy v^oi 

as i — >■ oo for the former and xi — ?■ for the latter. From the second equation of ([231 



X2 



yoi/ Woi Vol J 



(24) 



o Case 3: yo2 7^ 0. Similar to the case above, we get 

I X2 + ( 



xi = ( ) a;2 + ( — ""3" ) X2 + 0{xl) as X2 — )• 0. 



. 2^02 2^02 

For a more concrete example, consider the system 



xi \ /—I \ / Xi\ / xf + 8xiX2 + X2 



\X2 J 

Using ([24l) . provided yoi 7^ 0, 

'2/02 



2^2 



yoi 







Xl + 



t/ \x2l I 8x^ + X1X2 + 8x| 



yo2 
yol 



+ C'(xf) as Xl ^ 0. 



(25) 



(26) 



We can use the sign of (yo2/yoi)^ — 8 in ([26]) to deduce concavity. A phase portrait for ([25]) is 
sketched in Figure [TJ 



2.4.2 Example of Finding the Transformation 



Let 

and consider the system x = Ax + b(x). We will use the approximations (|2ip to find i/'(xo). Here, 



^(i)(xo)=xo and ^('"^^^(xo) 



xo2 + Ce^' [xr(«,V'(^'"Hxo)) 
An easy inductive argument and Proposition 1191 establish 



ds 



for 771 G N. 



(m+l)/ 



, a;o2 + x\ 



, and V(xo) 



X02 + 
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Figure 1: (a) The phase portrait near the origin for the example system (|25|) . (b) The sign of (2/02/yoi)'^ ~ 8 
in different regions of the j/oi2/o2-pla'rie, using (|15l) and (I26p . specifies concavity. 

2.5 Transforming Widely-Spaced to Closely-Spaced Eigenvalues 

This section deals with the case where the eigenvalues are closely-spaced relative to the nonlinear 
part of the differential equation. That is, k < a. Momentarily, we will deal with the case where 
the eigenvalues are widely-spaced relative to the nonlinear part. That is, k> a. However, the 
techniques for widely-spaced eigenvalues are decidedly more tedious. Fortunately, in practice it is 
sometimes possible to transform a system with widely-spaced eigenvalues to a system with closely- 
spaced eigenvalues. The transformation need not be polynomial or even analytic. 

The reduction of widely-spaced eigenvalues to closely-spaced eigenvalues is related to the theory 
of normal forms. Assume the nonlinear vector field b is analytic (or sufficiently smooth) at the 
origin and the eigenvalues of the matrix A are widely-spaced relative to b. 

Suppose that there is no resonance in the eigenvalues of A. For a given m S N with m > k, by 
the Normal Form Theorem there exists a sequence of near-identity, analytic transformations such 
that the differential equation x = Ax + b(x) is converted to a new differential equation of the form 
y = Ay + c(y) where ||c(y)|| = 0(||y|j™) as ||y|| — > 0. The new system, of course, has eigenvalues 
which are closely-spaced relative to the nonlinear part. 

Suppose now that there is resonance in the eigenvalues of A, say with lowest order m € N. The 
best we can do is find a sequence of near-identity, analytic transformations to convert x = Ax + b(x) 
to y = Ay + c(y) where ||c(y)|| = C)(||y|j'") as ||y|| — > 0. If m > k, then the eigenvalues of the new 
system will be closely-spaced relative to the nonlinear part. 
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3 Widely- Spaced Eigenvalues 
3.1 Introduction 

For simplicity, we will restrict our attention to a special case which is most applicable to us. The 
general case is outlined in §3.51 Consider the initial value problem 

i = Ax + b(x), x(0)=xo, where A:=(" ^) and b(x) = ( ) (27) 

\0 KaJ V^2(x)y 

with a < and k > 1. The matrix A and vector field b satisfy the same assumptions as before (see 
([3]) and (jH)). Additionally, we assume 

K>a. (28) 

That is, the eigenvalues are widely-spaced relative to the nonlinear part. Note that Ka < aa < a < 0. 
Note also that (f27|) can be written, where xq = (xoi,xo2)"^, as the integral equation 

cb(t^)=( + ^0 e-(*-^)6i(0(s,xo)) ds \ 

"P^ '^0^ l^e-»*a;o2 + /o e''«(*-^) 62(0(5, xq)) ds) ' ^ ' 

We can specialize the basic estimates to the case at hand. First, we know from Claim [3] that 
||e*^||=e"* and ||e-*^|| = e"''"* for all t > 0. (30) 
Furthermore, since A is diagonal, we know from ([6]) that there are S,k > such that 

||</)(t,xo)|| < A;e"* llxoll for all (t,xo)G05. (31) 
Claim 20. For any a > 0, there are 6,k > such that 

|(/)2(t,a^)| < A:e("''+'")* lla^ll for all {t,xo)en5- (32) 

Proof: Applying the estimates (HD and ([51]) to the second component of ([2^]) . there are 5,ki,k2 > 
such that if {t, xq) € il-s then 

I </'2 (t , xo ) I < fci e'^'^* 1 1 xo 1 1 + A;2 e^^' | e['^-(''-°)«]« ds | 1 1 xq 1 1" , 

where we inserted arbitrary o" > to allow for the possibility that k = a. Simplify to get (|32p . □ 
Remark 21. When k > q, we can take cr = in (j32p and throughout this section except for ^3.51 
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3.2 Iterates 

3.2.1 First Component of the Transformation 

In the previous section, where the eigenvalues were closely-spaced relative to the nonlinear part, 
we were able to "flip the integral" in the general integral equation ([2]). However, we cannot do this 
in this section. Why? The estimates ([4]), (pOl) . and (fSTl) tell us there are 6,k > such that 

||e-*^b(0(t,xo))|| < fce^"-")"* llxoir for all (t,xo) G ^s- 

The condition (I28p prevents us from concluding that ||e~*'^b(0(i, xo))|| decays exponentially. This 
in turn means that we cannot, in general, "flip the integral" in ([2]). Hence, we cannot find an 
appropriate yg and define the iterates as easily as with the case of closely-spaced eigenvalues. 

Claim 22. There are S,k > such that 

|e-'^*6i(0(t,a^))| < A;e("-i)"* ||a;oir for all it,Xo)ens. (33) 

Proof: It follows from (g]) and i^. □ 
Since a > 1, we see from (I33p that we can "flip the integral" in the first component of (129p . 

Claim 23. There is a 5 > such that 



e 



at 



M^)-ir e''^'-'^bi{(l,{s,xo)) ds^ 



cf>it,Xo)=i ^''^' -^y^^-^-^^^ - \ (34) 

\e^'''xo2 + foe^''^'-'^b2{^{s,xo))ds J 

for all {t, Xo) € ^Is, where 

iji{xo) := xoi + e-''%i{cf){s,xo)) ds. (35) 
Jo 

Claim 24. There are 5,k > such that 

\iIji{xo) — xoi\ < k \\xo\\"' for all xo ^ Bs- (36) 

Proof: It follows from ([33]) and ([35]). □ 
3.2.2 First Group of Iterates 

The iterates, which will be denoted by {^'•'"^(t, yo)}^^]^) must be split into two groups. Define 
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which is such that the p*^, but not (p — 1)* , iterate will be close enough to (f){t,xo) so that we can 
perform another "trick" in the second component of ()34p . It can be shown 



pe{2,3, ...} and a + {p - 2) /3 < k < a + {p - I) f3. (37) 
If llyoll is sufficiently small, we take as the first iterate 



and (if p > 2) define recursively 



To connect these iterates with the flow, choose 7/01 '■= ''/'i(xo) with ||xo|| sufficiently small. 

Remark 25. Since the first group of iterates does not depend on yo2) we write x*-™^ (t, yoi) instead of 
X^^\tiyQ)- This is necessary since V'2(xo) will be defined in terms of X^^~'^\t^ V'i(xo))- However, if 
we need to refer to all the iterates at once we will use, for simplicity, the notation {x*-™^(t, yo)}^_^- 

3.2.3 Existence, Decay Rate, and Closeness to the Flow of the First Group 

Proposition 26. Consider the transformation defined by (j35p and the first group of iterates defined 
by ()38p . Then, there is a 6 > (independent of m) such that if \\xq\\ < 6 then ^^"^^(t, 2/01) exists 
for each m € {1, . . . ,p — 1}, where yoi := ^i{qcq). Moreover, there is a k > (independent of m) 
such that 

X^"'Ht,yoi)\\<ke'''\\xo\\ for all {t, xq) e Qs, m e {1, . . . ,p - 1} . (39) 
Proof: The proof is similar to that of Proposition [T3] and omitted in the interest of space. □ 

Theorem 27. Consider the transformation defined by (|35p and the first group of iterates defined 
by ([38]) . Let a > be sufficiently small. Then, there exist 6 > and {km}^=i C such that 

0(t,ato) -X^™^(t,yoi)|| < A:™e[-+(™-i)'^l["+-]*|l^ll (40) 
for all {t, qcq) € $^5 and m S {1, . . . ,p — 1}, where yoi '■= V'i(^)- 

Proof: The proof is similar to that of Theorem [15] and omitted in the interest of space. □ 

Remark 28. Notice in (|40p we have |jxo|j as opposed to ||xo|j"^'-™ "'^^'^ which we had in (|20p . This 
is the result of e'^"*2;o2 in ()34p having no corresponding term in ()38p with which to cancel. Note 
that we cannot merely insert e'^'^*xo2 into the definition of the second component of the iterates 
since we need the iterates to be defined independently of the initial condition xq. 
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3.2.4 Second Component of the Transformation 



Let uqi := ^i(xo), where ||xo|| is sufficiently small. Consider the integral equation (|34p for which the 
first integral has been "flipped" but the second integral has not. We will manipulate the expression 
for the second component of the flow so that we can "flip some integral" in some way. Write 

ft 



^Kat 



+ e' 



Kat 



X02 + 
t 



ds 



(41) 







[62(0(s,xo))-62(x(P-i)(s,yoi) 
''b2(x^^-^\s,ym)) ds. 

Before we "flip the first integral" in (|41|) . we need to check the decay rate of the integrand. 
Claim 29. Let a > be sufficiently small. There are 6,k > such that 

e-""* [b2m,xo)) - b2{x^P-'Ht,yoi))] \ < k e^"+(P-^')''-->+^^' \\xo\f+' (42) 
for all {t, Xq) € where yoi := ipi{^)- 

Proof: It follows from (0]), pip . (j39p . and (j40p in conjunction with the Mean Value Theorem. □ 
Claim 30. There is a 6 > such that 



for all {t, xq) G Qs, where 

i^2{Xo) ■= X02 + 



e-'^'M^)+ f e'^'^(*-^)62(x(^^')(5,V'i(^))) 



ds 



Ka(t—s) 



b2{4>{s,Xo)) - b2(x''P-^Hs,Mxo)) 



b2{cl>{s,Xo)) - b2(x^P'^Ks,M^)) 



ds 



ds. 



Proof: It follows from (|34l). (l37l). (j4T|). (|42|). and 

Claim 31. There are 5,k > such that 

|V'2(a^) - 2;o2| < A; ||a^||''+"^ for all xq e Bs- 

Proof: It fohows from (gS]) and (gH). 

Notice the definition of ^^2(^0) depends on ^i(xo). See ^3.5.21 for how to approximate 

xoi + /o°° e""''6i(0(s,xo)) \ 



V'(xo) 



X02 + /c 







[62(0(5, xo)) -62U(^-'Hs,V'i(xo)))] ds 



(43) 

(44) 
□ 

(45) 
□ 

(46) 



Claim 32. There are 5,k > such that 



\\tl}{xo) — xqW < k\\xQ\\"^'"^"'^^^^ for all xq e Bs- 
Furthermore, xj) is a near-identity transformation. 

Proof: The conclusion follows from (j36p and (j45p along with the fact that a > 1 and /3 > 0. □ 
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3.2.5 Second Group of Iterates 



We are now in a position to define the remainder of the iterates, {x*-™'^^, yo)}^_p- Let ||yol| be 
sufficiently small. With inspiration from (1340 and ()43p . we will define the p^^ iterate by 



2"*yoi e^^'-'^bi{x^P-'\s,yoi)) ds' 



and define the subsequent iterates recursively by 



(47a) 



{p)l 



b(x('")(s,yo)) -b(x(P-i)(.,yoi))] ds (m > p) . 

(47b) 

To connect the iterates with the flow, we choose yg := i/'(xo), where ■0 is defined in (|46p and ||xo|| 
is sufficiently small. Observe that we can combine (j34p . (j43p . and (j47ap to obtain the following. 

Claim 33. There is a 5 > such that, for all (t, Xq) G il^, 



^{t,xo) = x'^P\t,'iP{xo)) 



bicf>{s,Xo))-b(x^P'^Ks,M^)) 



ds. 



3.2.6 Existence, Decay Rate, and Closeness to the Flow of the Second Group 

Proposition 34. Consider the transformation defined by ()46p and the second group of iterates 
defined by ([I7|) . Then, there is a 6 >0 (independent of m) such that if \\xo\\ < 6 then x^"^Ht,il}{xo)) 
exists for each m > p. Moreover, there is a k > (independent of m) such that 

%(")(*, -0(a3D)) <A;e"*||a^|| for all {t, xq) e n^, m > p. 
Proof: The proof is similar to that of Proposition [T3] and omitted in the interest of space. □ 

Theorem 35. Consider the transformation defined by ([I6|) and the second group of iterates defined 
by (|47p . Let a > be sufficiently small. Then, there are constants 6 > and {km}m=p ^+ -^^c/i 
that, for all (t, Xq) € and m > p. 



ct>{t,x,)-x^^\tM^)) 



Proof: The proof is similar to that of Theorem [T5] and omitted in the interest of space. 



□ 



Corollary 36. Consider the iterates defined by (|38p and (|47p . There exists 6 > such that 
lim x^™''(t, ■0(2^)) = 0(t,a^) uniformly in Qs- 



Proof: The proof is similar to that of Corollary [T71 



□ 
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3.3 Resonance, Quadratic Nonlinearity, and Relating Components 

The special case which we focus on here is applicable to the Michaelis-Menten mechanism which 
we will see later. Assume k G {2,3, . . .}, that is, there is resonance in the eigenvalues. Define 



b(x) := 



binxf + 6112x1x2 + 6i22x| 
6211X1 + 6212x1x2 + 6222X2 , 



where {bijk}f j C M are constants with at least one bijk being non-zero. Take a = 2, /3 = 1, and 
p = K. Hence, consider the initial value problem 

X = Ax + b(x) , x(0) = XQ. (48) 

Let Yo := '0(xo), where ||xo|| is sufficiently small, satisfy yoi > 0. (Since xi{t) = yoi + o(e'**) as 
t 00, we assume yoi > since it corresponds to solutions with xi > approaching the equilib- 
rium point in the slow direction. This is desirable since it corresponds to, for example, solutions 
that arise in chemical kinetics.) Also, let x(t) be the solution to (|48p . Our goal is to obtain an 
asymptotic expression for X2 in terms of xi as xi — > 0"*" up to an order for which the initial condition 
distinguishes solutions. This can be achieved without knowing yg in terms of xq explicitly. 

o Case 1: k = 2: We will only need the first two iterates. Using (|38ap and (|47ap . 



V y U2e^'^*+[62iiyo^ 



g2at 



Claim 37. The components xi{t) and X2{t) of the solution to (j48p satisfy 
xi(t) = (yoi)e'^* + o(e'^*) and X2(0 = (62iiyo'i) * e^'^* + (^02) e^"* + o(e2'^*) as t ^ 00. (49) 

Proof: It follows from Theorem [35] and our expression for the second iterate. □ 
Claim 38. We can write e'^*, e^°*, and t in terms of xi when xi 0^ as 

e-*= f— ^xi + o(xi), e^-i = (^)xl + o(xl), and t =(-) lnixi)-i^^^ + 0(1) . (50) 
\yoiJ KVoiJ \aj a 



Proof: To obtain the first part, use the first expression of (|49p and a method employed in 
^2.4.11 The second part follows from the first. To obtain the third part, write the first part 
as e"^* = [1/2/01] [xi] [1 + 0(1)] as xi — )• 0"*". Taking the natural logarithm of both sides, 

at = -ln(yoi) + In(xi) + ln(l + 0(1)) as xi — 0+. 

Observing that ln(l + u) = o(l) as u ^ 0, we are left with what we were trying to show. □ 
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Proposition 39. Consider the initial value problem (|48p . where k = 2, \\xq\\ is sufficiently small, 
Vo ~ c-''^^ V'i(^) > 0. If x{t) is the solution, then we can write X2 in terms of xi as 



X2 



'J211 



x\ ln{xi) + 



yo2 b2nln{yoi) 



xl + o{xl) 



as xi 



Proof: Substitute the second and third of (|50|) into the second of 



0+. 

then simplify. 



□ 



o Case 2: k£ {3, 4, . . .}: Using (pHD . 



. 



at 



and X^^\t,yoi] 



yoi e"^ + 



a 



J2at 





g2at _ 


^2 11^01 


(2-K)a 




(2-K)a 



^Kat 



By Theorem 1271 



xi(t) = yoie"*+('^iii^)e2'^* + 0(e3'^*) and x^it) 



(2 - K)a 



3at\ 



as t — 7- oo. By inspection, X2 = 0(^xfj as xi ^0^. To take this further, we need the following. 

Claim 40. There exist constants and {pi}'^Z2 (which do not depend on the initial condition) 

and a constant g (which may depend on the initial condition) such that, as t ^ oo, 

K 

Xl{t) = Y,^^{e'''yQlY + o{e^^^^^-'), with 6 = 1, 

h2ii 



i=l 

K-1 



and X2(t) = Xl/'i(e"*yoi)' + £'(e""*)+o(e("+^^"*), with p2 



(2 - K)a' 



(51a) 
(51b) 



Proof: It follows from Theorems [27] and [35} the definition ()38p of the first k — 1 iterates, and the 
definition (j47ap of the k)'^ iterate. □ 

Proposition 41. Consider the initial value problem (148 p . where k € {3,4, . . .}, ||a:^|| is sufficiently 
small, Uq = iP{xq), and i/;i(ato) > 0. If x(t) is the solution, then for some constants {ci]1Z2 (which 
depend only on the differential equation) and C (which depends on the differential equation and the 
initial condition) we can write X2 in terms of Xi as 



X2 = ^^Cix\ + C xl + 0[x1'^^) as xi^Qi^, with C2 



b 



211 



(2 - k) a 



(52) 



Proof: Note that xi = ©(e***) as t — > oo and e"^* = 0{xi) as xi ^0"*", which follow from ()51ap . 
Now, by inverting (j51ap in a process we have done many times before, it follows that there are 
constants {i^ilJLi 'which do not depend on the initial condition such that 



e'^*yoi = ^' ^1 + C'(x5'+i) as xi ^ 0+, with z^i = 1. 



(53) 



i=l 
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Substituting (l53|) in (l5Tbj) . we are left with (j52p for some constants {ci}^^2 (which do not depend 



on the initial condition) and C (which may depend on the initial condition). 



□ 



3.4 Generalizing to a Diagonalizable Matrix 

Suppose A and b are as before except A is only assumed to be diagonalizable with (fast and slow, 
respectively) eigenvalues Ai = a and A2 = na. Consider the initial value problem 



X = Ax + b(x) , x(0) = xq. 



(54) 



Now, there exists invertible matrix P such that A = P ^AP, where A := diag (a, ko). Define 
r(u) := P~^b(Pu) and uq := P^^'^xq. The following is easy to verify. 

Claim 42. If x{t) is a solution of (I54p . then u{t) := P^^x{t) is a solution of 

u=Au+r{u), tx(0) = ■Uq. (55) 
Moreover, there are 6,ki,k2 > such that ||r(M)|| < ki and \\Dr{u)\\ < k2 \\u\\^ for all u € Bs- 



3.4.1 Associated Iterates 

The iterates {x^™^(t, ■'/'(uo))}^^^^, defined with respect to A and r, allow us to construct associated 
iterates, {Px^"^) (t, ■0(P~^xo)) }^^-^, for the original system. 

Theorem 43. Let x{t) be the solution of the initial value problem (I54p . Consider the iterates 
given by (j38p and ()47p and the transformation given by ()46p . defined with respect to ()55p . For any 
sufficiently small o" > 0, there exist 5 > and {km}m=i ^+ such that, for all {t, qcq) € il-s and 
m G N, 



x{t)-Px^"'\t,'4,{p-'xo)) 



^ I, 11^ |,max{0,-m+l-p}/3+l 

\ Km. S 1 1 Xq 1 1 



Proof: Let 5 > be sufficiently small so that if ||xo|| < 5 then the appropriate estimates apply. 
Hence, assume xq € Bs- Let u(t) := P^^x(t), which is the solution to (|55p . Observe that 



uqII < IIP """11 ||xo|| and 



x(t)-Px('")(i,./,(uo)) 



< IIPI 



u(t)-x('")(t,^(uo)) 



Applying Theorems [27] and [35] give the conclusion. 



□ 



3.4.2 Resonance, Quadratic Nonlinearity, and Relating Components 

Consider the initial value problem (|54p with the particular nonlinear part 



b(x) : = 



biuxj + 6112x1x2 + bi22xl 
b2uxj + 6212x1x2 + 6222^1 



(56) 
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where {bijk}" 



C M are constants with at least one h 



I i,j,k=l 

the diagonalized initial value problem ()55p is of the form, where {''ijfcli ^ ,^=1 constants, 



•ijk being non-zero. The nonlinear part for 

2 



r(u) 



nuuj + ru2UiU2 + ri22ul 
r2iinf + r2i2UiU2 + r222u\ 



We want to be able to write X2 as a function of x\ when x(t) approaches the origin in the slow 
direction and there is resonance in the eigenvalues, that is, when k G {2, 3, . . .}. (The no-resonance 
case K ^ {2, 3, . . .} can be handled using iterates but also by using Poincare's Theorem. See, for 
example. Theorem 19 of the authors' [6].) Observe that the material of ^3.31 applies to (f55|) . 

Since x(t) approaches the origin in the slow direction from the right, we can assume pn > 
and ^01 > 0, where vq := '0(uo) with ■0 being defined with respect to the diagonalized system (f55|) . 
To see why we can take pn > 0, observe that we can take the first column of P to be the slow 
eigenvector and we can take the slow eigenvector to have positive first component. Note the slope 
of the slow eigenvector is "p^ilvwi so we expect X2 = {P2i/pii) xi + o{xi) as xi 0+. To see why 
we can take vqi > 0, observe 

ui{t) = Vol e"^ + o(e''*) , U2{t) = o(e''*) , and xi{t) = pnvoi e"^ + o(e''*) 

as t ^ oo, the last of which uses the relationship xi = puui + pi2U2- 

Theorem 44. Consider the initial value problems (I54p . with b{x) given by ()56p . and (I55p . Suppose 
K £ {2,3, . . .}, pii > 0, \\xq\\ is sufficiently small, and vqi > where Vq := t/^iuo) with ijj as in (jl6|) 
defined with respect to ([55]) . Let x{t) be the solution to (j5^ . 

(a) If K = 2, then we can write X2 in terms of xi as 

'vq2 r2iiln{piiVQi) 



X2 = 


'P21 


Xl + 


V2II \P\' 


xl ln{xi) + 


'\P\' 






.Pll. 




. Pna . 


.Pii. 





^01 



xl + o(xl) 



(b) /f K G {3, 4, . . then we can write X2 in terms of xi as 

xi + ^Cix\ + Cx1 + 0{x1'^'^) as xi-j-O^ 



X2 



P21 



Pii 



as Xl ^ . 

(57) 



(58) 



i=2 



for some constants {ci}'^^2 (which depend only on the differential equation) and C (which 
depends on the differential equation and the initial condition). 



Proof: As mentioned above, x(t) approaches the origin in the slow direction and is strictly positive 
for sufficiently large t. 
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(a) Assume we can write U2 as a function of ui, say U2 = f{ui). Then, xi = puui + pi2f{ui) since 
X = Pu. Assume further that we can solve this relation for ui in terms of xi, say ui = g{xi). 
Thus, since x = Pu we have X2 = h{xi), where h{xi) := P2ig{xi) + P22fig{xi)). With this in 
mind, the relation ()57p . after routine simplification, follows from Proposition 1391 and techniques 
already employed in this paper. Note that |P| = puP22 — P12P21 7^ 0. 

(b) The relation (|58p follows from Claim HOl the fact x = Pu, and the method used in the proof 
of Proposition [5TJ 



Remark 45. Suppose that, in Theorem the components of b(x) are polynomial (instead of 
purely quadratic) with lowest order a € {2, 3, . . .}. If k = a, then the expression for X2 in terms of 
xi as xi — >■ 0+ will contain In(xi). If k > a, then the expression for X2 will be a Taylor series in xi 
with Cxi being the lowest-order term which depends on the initial condition. 

3.5 Generalizing to Higher Dimensions 
3.5.1 Introduction 

The techniques of this section can be generalized to the n-dimensional case where the eigenvalues 
of A are widely-spaced relative to b(x). The derivations and proofs will be omitted since they 
are tedious and reminiscent. Consider the general initial value problem ([1]) with n > 2 and k > a. 
Suppose that J G M"^" is the real Jordan canonical form of A which satisfies the same ordering ([3]) 
for the eigenvalues {Aj}^^^. Then, there exists invertible matrix P G M"-^"- such that J = P~^AP. 

We need to partition the Jordan blocks of J. Suppose the eigenvalues of J have ^ € {2, . . . , n} 
distinct real parts {p-jYj^i, each with multiplicity nj so that X]j=i "-j = which are ordered 
according to /i^ < • • • < /Ii < 0. Then, we can write J = ©j=i Jj where each Jj G is the 

direct sum of all real Jordan blocks having associated eigenvalues with real part Jlj. Note that 
Jli = /Ui and fli = fin- We will write all vectors x G in component form as x = (xi, . . . ,X£)^ 
with Xj G for each j G {1, . . . , i}. 

Hence, consider the initial value problem 



where u := P x, uq := P xq, and r(u) := P b(Pu). Denote the flow of ([^^ by ip{t, uq) and 
observe each component satisfies the integral equation 



□ 



u = Ju + r(u) 



u(0) = uo 



(59) 




Define 



Kj : — 



^ (j = l,---,^) and jo := m\n{ j e {1,... ,i} : Kj > a} 
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Component j 
Index m 


Definition or x) (^jVq) 


1 < J < jo 
m — 1 




1 < j < jo 
TO > 2 




jo <j<i 







TO = 1 






jo <j<e 

2 < m < pj 




Jo <j<i 

TO = Pj 




JO <j<i 

TO > Pj 




ds 



Table 1: Definition of tlie iterates {x'^™''(^; '^o)}^_]^ for tlic initial value problem (|59p . 



and note 1 = /ti < K2 < • • • < = k and jo € {2, ... , £}. Observe 1 < j < jo corresponds to blocks 
which are closely-spaced relative to the nonlinear part. Define also 



Pj 



1, if 1 < j < jo 

+2, if jo<J< 



(j 



Observe 2 < pj^ 



< 



< p£ and 

a + {pj -2)f3<Kj <a + {pj - 1) /3 (j = jo, 



Note that it can be easily shown, for any o" > 0, there exist 5 > and {kj}.^.^ C such that 
uo)|| < kj e(°^'i+°')* ||uo|| for all (t, uo) G 0^ and j G {jo, . . . ,^}. This estimate is required to 
establish the base case of Theorem 1491 below. 



3.5.2 Iterates and the Transformation 



Assume that ||vo|| is sufficiently small so that we can define the iterates {x^^^i, vo)}^_-^. The 
definition, given in Table [H must be split into many cases. 

Remark 46. For a particular j G {jo, . . . ,^}, voj only appears in the iterates {x^^H^' '^o)}^ 
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Component j 
Index m 


Definition of ■0j™^(uo) 


1 < J < ^ 
m = 1 




l<j < jo 
m > 2 




) ds 




Jo <J<i 
m > 2 





(s,tA(™-i)(uo)))' 


ds 



Table 2: Definition of tlie approximations |'0'^™^(uo)}^_-|^ to tlie transformation '4>{uq) for tlie initial value 
problem (ICTl) . 



To connect the iterates with the flow, where ||uo|| is sufficiently small, take vq := i/'(uo), where 

^ \ uoj + e-'^^r^{if{s, uo)) ds, if I < j < jo 
tpjiuo) := < • (60) 

(uo, + e-^-J^ [r,((^(5,uo)) - r, (x^^^-^H^, ^(uo)))] ds, if Jo <j<i 

It can be easily shown that, for some constants 6,k > 0, the transformation tp satisfies 

fiin{a,/3+l} 



||'0(uo) - Uoll < k ||uo| 



for all uq G iS^. 



Remark 47. The term x ('S)'0(uo)) inside the integral in (|60p for the case j'q < j < ^ does 
not depend on the components {i/'i(uo)}^^j. Thus, the definition is explicit. 

The iterates given in Table [H which are expressed in terms of the unknown parameter vq, 
are useful for obtaining the form of the asymptotic expansion for the flow (f{t, uq) as t — > cxd when 
Vq = i/'(uo). Fortunately, just like with closely-spaced eigenvalues, we can approximate ip{uo). The 
definition of the approximations {■0'''")(uo)}^_j^ is given in Table [21 



3.5.3 Results 

Proposition 48. Consider the iterates defined in Table{l\ with respect to the initial value problem 
()59p and let a > be sufficiently small. Then, there is a 6 > ( independent of m) such that if 
\\Uq\\ < 6 then 1/7(1*0)) exists for each m G N, where tp{uQ) is as in ()60p . Moreover, there is 

a k > ( independent of m) such that 



for all {t, Uq) G Qs, rn eN. 



(61) 



Theorem 49. Consider the iterates defined in Table{l\with respect to the initial value problem (|59f 
and let cr > be sufficiently small. Then, there are constants 5 > and {A;m}m=i ^+ ■'^^c/i that 



ct>{t,xo)-Px^^\t,'4,{p-'xo)) 



^ km S 



a+im-l)|3][^ll+cT]t ||^||max{0,-m+l-p4/3+l ^^^^ 
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for all {t, xq) G o-nd m £ N, where cj>{t, Xq) is the solution of the initial value problem ([T]) and 
iI^{uq) is as in (f60l) . 

Remark 50. If A is diagonalizable and > a then we can take cr = in (j6ip and (|62p . 

Proposition 51. Consider the transformation defined by (j60p and the approximations defined in 
Table\^ There are constants 5 > and {fcm}m=i ^ ^+ •5'uc/i that 



4 Application to the Michaelis-Menten Mechanism 

The Michaelis-Menten mechanism of an enzyme reaction will be our main application. Particularly, 
we will specify the asymptotics of the scalar reduction for small substrate concentration. 

4.1 Introduction 

The Michaelis-Menten mechanism, named after two researchers who worked on the model, is the 
simplest chemical network modeling the formation of a product from a substrate with the aid of an 
enzyme. See, for example, [101 HSl EZl IIHI EO] . There are different and more complicated networks 
involving, for example, multiple enzymes, inhibition, and cooperativity. However, even though the 
Michaelis-Menten mechanism is relatively simple, there are many results on the mechanism and the 
standard techniques of analyzing the reaction are easily adapted to the more complicated networks. 

In the scheme of Henri, later explored by Michaelis and Menten, an enzyme E reacts with the 
substrate S and reversibly forms an intermediate complex C, which decays into the product P and 
original enzyme. With reaction-rate constants A;_i, ki, and k2, this can be written symbolically as 



Technically, the reaction from C to P + E \s reversible, but the rate constant k^2 is usually negli- 
gible and the reverse reaction is omitted. Henri's original model allowed for the reverse reaction. 
Moreover, in a cell or laboratory the product is harvested preventing the reverse reaction. 

Let lower-case letters stand for concentrations and r for time. The Law of Mass Action gives 

ds de dc dp 

— = A;_ic — Kise, — = + K2) c — kise, — = kise — (/c_i + K2) c, — = k2C. (63) 

ar ar ar ar 

Here, the initial conditions are arbitrary. Traditionally, s(0) = sq, e(0) = eo, and c(0) = = p{0). 

Observe that, using (f63]l . the concentrations s(r), c(r), and e(r) do not depend on the value of 
p{t). Moreover, it follows from ([63]) that c(r) and e(r) satisfy the conservation law c(r) + e(r) = eg. 
We can hence consider only the system of two ordinary differential equations 



m/n{a,/3+l}+(m-l)/3 



for all iio G Bs- 



S + E^C%P + E. 



fe_i 




kis (eo - c) - {k-i + k2) c. 



(64) 
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This is frequently referred to as the planar reduction of the Michaelis-Menten mechanism. The 
traditional initial conditions are s(0) = sq and c(0) = 0. However, we will again allow the initial 
conditions to be arbitrary. The planar reduction (j64p can be reduced as well to the scalar reduction 



dc kis (eo — c) — {k_i + ^2) c 
ds k_ic — kis (cq — c) 

We will make one last transformation. Define 

kis c kieo - k2 

t := KiCqt, X := — , y := — , e := — , and 77 := 



(65) 



+ k2 eo k^i + k2 fc-i + k2 

which are all dimensionless, and note e > and < rj < 1. (We make no assumption on the size of 
e, which is traditionally assumed to be small.) Thus, t, x, and y are, respectively, a scaled time, 
substrate concentration, and complex concentration. Observe x = s/Km and e = e^/Km, where 
:= (A;_i + k2)/ki is the Michaelis constant (the half-saturation concentration of the substrate). 
It is easy to verify that the planar differential equation (j64p assumes the dimensionless form 

X = —X + {1 — r] + x)y, y = [x — {1 + x) y] , (66) 

where ' = d/dt. Moreover, we can also write the scalar differential equation (j65p in the form 

y' = x-{l + x)y ^^^^^ ' .= A (67) 

e [—X + (1 — + x) y] ' dx 

4.2 Behaviour at the Origin 

The planar reduction (j66p . along with initial conditions x(0) = xq and y(0) = yo, can be written 

X = Ax + b(x) , x(0) = Xq, where A := ( ^ \ and b(x) := xy i ) • 

The matrix A has two distinct, real, negative eigenvalues given by 



-{e + l)±^{e + lf -Aeri 



2e 

Here, A^ is the slow eigenvalue and A_ is the fast eigenvalue. It can be shown that the eigenvalues 
satisfy A_ < — 1 < —r] < A_^ < 0. Moreover, it can be shown that the nonlinear part b(x) satisfies 



b(x)|| < ( Wl + e-2] ||x||2 and ||Db(x)|| < ( VT+T^) ||x 



Define 



A, +1 , A e + l + Jie + lf-4eij 
and K := — = . 

^ " \ e + 1- J{e + if - Aerj 
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Here, and a_ are, respectively, the slopes of the slow and fast eigenvectors of A. It can be 
shown a_ < 0, 1 < < 1/(1 — r/), and k > max [e, > 1- 
We want to apply Theorem 1441 Consider the matrices 



' 1 1 



a 



and A 



P-^AP 





A 



If we define u := P ^x, uq := P ^xq, and r(u) := P ^b(Pu), then the initial value problem ([68 
is transformed into an initial value problem of the form considered in ^3.4.21 Note that 

^ ^ cj, n? + (cr, + (7_) M1M2 + cr^Uo / -(T_ - \ ^ o" , (cr , + 

r(u) = ^ I I and r2u - ' + 



where r2ii is the coefficient of the uf term in the second component of r(u). 

The traditional approach of finding the behaviour of a scalar solution y{x) to ()67p as x — )■ 0^ 
is to obtain a power series solution. Using the ansatz y{x) = X^^q"^"^" ([67|) yields 

0, if n = 



(7^, if n = 1 

^n-l 



T.k=2 [in-k)a„-k + i'^-v) (n-fc+l)o-„,fc+i](jfc+[(n-l)q-i+g ^] an-i 

e~^ + (l— ry) (n + 1) (7i— n 



if n > 2 



(69) 

One problem with using the resulting series is that the denominator in the recursive expression for 
cr„ in (j69p is zero if and only if n = k. The authors showed this in [6j. This issue is resolved in the 
following theorem, which generalizes Lemma 22 of the same paper. 

Theorem 52. Let x{t) be a solution to the initial value problem (I68p . where \\xq\\ is sufficiently 
small, which approaches the origin in the slow direction from the right. Lety{x) be the corresponding 
scalar solution to (j67p and consider the coefficients {(T„}^-|^ given in (jGOp . Then, we can write 

\_k\ 



y{x) 



^a„x" + Cx'^ + o(x''), if K0{2,3,...} 



n=l 

a^x + 

K-1 



x"^ ln{x) + C + o(x^) , if K 



as X 



0+, 



^(T„x" + Cx'^ + 0(x'"+^), if KG {3,4,...} 



^ n=l 



where C is some constant (which depends on the differential equation and the initial condition). 



Proof: The first part was addressed in Lemma 22 of the aforementioned paper. Note the coeffi- 

I K 1 

cients of the integer powers of x must indeed be {an}!^^^ since they are generated uniquely by the 
differential equation. The second and third parts are immediate consequences of Theorem HH □ 
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The authors showed, in [6], that for any e > there is a unique attracting solution y = M.(x) 
(the slow manifold) to (I67p which satisfies 



H{x) < Ai{x) < a{x) for all x > 0, where H{x) :- 



X 



1 + x 



and a{x) :- 



X 



are, respectively, the horizontal isocline (the quasi-steady-state manifold) and the isocline for slope 
. Theorem [52] tells us the asymptotic behaviour of the slow manifold (and other solutions) at 
the origin up to the term for which solutions are distinguished by initial conditions. 

Planar solutions approach the slow manifold exponentially fast and then follow it as time pro- 
ceeds. Hence, it is of value to study the dynamics on the slow manifold which is of lesser dimension 
than the original system. Moreover, knowing the position of the slow manifold accurately allows, 
for example, all rate constants to be determined independently from steady-state data |23j . 

Remark 53. We can use Theorem [52] to comment on the validity of two hyperbolic rate laws, 
specifically the quasi-steady-state approximation H(x) and the a-approximation a{x). For any 
physically-relevant initial condition xq (that is, xo,yo ^ and xq 7^ 0) to the initial value problem 
(j68p . as shown in [6], the solution x(t) satisfies x(t) ^ as t — > 00 with the approach being in the 
slow direction from the right. Since H{x) = x + O(x^) and a{x) = a^x + O(x^) as x — )• 0^ and 
x{t) = 0^e^+*^ as t ^ 00, it follows from Theorem [52] that, in terms of the dimensional quantities 
of the original planar reduction ([M]) . if so,co > and (so,co) 7^ (0,0) then 



c(t) 



eos(T) 

Km + s{t) 



+ e 



fcieoA , 



eos(r) 



Km(T + s{t) 



+ < 



'0(^e['^i^"^-^) , if 1<K<2 

0j^^g[2fcieoA^]r^ ^ if K = 2 

o(e[2'=i"°V]") , if K>2 



as r —7- 00. In all cases of the size of k, the a-approximation is better than the quasi-steady-state 
approximation. (In fact, there is no better hyperbolic rate law than the a-approximation.) 



5 Open Questions 

It would be useful to allow the origin to be degenerate, that is, replace /xi < in ([3]) with /xi < 0. 
This change would not be trivial because, generally, the origin would not be asymptotically stable 
(there may be hyperbolic sectors) and there would not be exponential decay for the flow. Moreover, 
it would be an interesting exercise to specify the asymptotic expansion of solutions for a specific 
system when the initial condition is on an invariant manifold. Finally, similar to Remark [53] 
the rate of exponential-time decay for the error in the semi-infinite time interval version of the 
Quasistatic-State Approximation Theorem in §8.3 of [16] should be determinable. 
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